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ABSTRACT 

Mlnisonde  data  collected  in  the  Athabasca  Oil  Sands  area  from 
1975  to  1979  were  analysed  to  determine  regional  values  of  roughness 
length  (Zq).  A rigorous  selection  procedure  reduced  the  working  data 
set  to  a small  fraction  of  the  original  size.  A least  squares 
technique  was  used  to  determine  Zq  from  profiles  of  wind  and 
temperature  typically  measured  near  the  50,  100,  and  150  m levels. 

Mean  Zq  values  calculated  with  allowance  for  diabatic  and 
displacement  height  effects  ranged  from  about  8 m downwind  of  the 
Syncrude  plant  site  to  about  1 m in  the  Athabasca  River  valley. 
Uncertainties  in  the  estimates  were  of  the  same  magnitude  as  the  mean 
values.  No  differences  in  Zq  were  found  with  wind  direction.  The 
large  values  for  Zq  were  attributed  primarily  to  form  drag  from 
terrain  features  In  the  area  during  slightly  unstable  conditions.  The 
study  suggested  that,  where  form  drag  Is  important,  Zq  may  be 
stability  dependent. 

An  error  analysis  using  reasonable  uncertainties  for  wind 
speed,  balloon  height,  and  temperature  gradient  measurements  showed 
that  probable  errors  in  the  estimate  of  Zq  were  comparable  to  the 
observed  variability  in  Zq. 
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INTRODUCTION 


1 . 


The  roughness  length  is  a parameter  relating  wind  drag  at  a 
surface  to  the  vertical  gradient  of  horizontal  wind  speed.  It 
indicates  the  extent  to  which  mechan i ca I I y- i nduced  turbulence  is 
generated  by  wind  flow  over  particular  surfaces.  For  air  quality 
studies,  the  roughness  length  enters  the  dispersion  formulation  either 
directly  or  by  means  of  site-specific  empirical  parameters.  For 
example,  in  the  Gaussian  frequency  distribution  model  developed  by 
Davison  et  al.  (1981a),  the  roughness  length  (Zq)  is  used  explicitly 
to  calculate  the  friction  velocity  and,  hence,  the  fluctuations  of  the 
wind  for  the  dispersion  formulation. 

The  general  study  of  wind  profiles  and  surface  stress,  of 
which  measurement  of  roughness  length  is  a part,  has  proceeded  from 
uniform,  flat  terrain  to  various  types  of  topography  and  surface 
conditions.  Walmsiey  et  al.  (1982)  and  Taylor  et.  al.  (1983)  described 
the  application  of  a three-dimensional  numerical  model  based  upon  Mason 
and  Sykes  (1979)  to  a small-scale  terrain  feature  (Kettles  Hill, 
Alberta).  An  extensive  low  level  wind  profile  measurement  program  at 
Kettles  Hill  (Taylor  et  al.  1982)  provided  encouraging  experimental 
support  for  the  model  predictions.  However,  the  results  indicated  a 
need  for  spatially  varying  roughness  lengths. 

Generally,  the  effective  roughness  length  has  been  found  to 
be  both  direction-  and  height-  dependent.  Beljaars  et  al.  (1983)  found 
that  the  friction  velocity  (u*)  changed  with  height  in  response  to 
changes  in  upstream  roughness  and  that  the  vertical  velocity 
fluctuations  tended  to  scale  with  a "local”  u*  at  3.5  m whereas 
horizontal  fluctuations  tended  to  scale  with  global  values  of  u* 
measured  at  22.5  m.  Ming  et  al.  (1983)  analysed  routine  wind  profiles 
from  three  towers  over  100  m h i gh  in  various  types  of  complex  terrain 
in  New  England.  They  calculated  effective  roughness  lengths  for  upper 
and  lower  parts  of  the  profile  and  attempted  to  relate  the  values 
obtained  to  surface  features  various  distances  upstream  using 
Hojstrup’s  (1981)  relationship  between  height  on  a tower  and  upstream 
distance  of  influence.  Ming  et  al.  attributed  some  of  the  large 
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roughness  lengths  in  the  upper  layers  (as  large  as  11  m)  to  form  drag 
effects  due  to  low  hills  (100  to  200  m)  upstream  of  the  tower  site. 

In  the  present  study,  minisonde  data  from  the  Athabasca  Oil 
Sands  area  were  used  to  estimate  a value  of  roughness  length 
appropriate  to  the  region.  The  minisonde  profiles  were  collected  and 
processed  by  various  groups  and  made  available  by  Alberta  Environment 
as  a digital  data  set. 
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2.  IMPORTANCE  OF  THE  ROUGHNESS  LENGTH  IN  THE  ATHABASCA  OIL  SANDS 


Mechanical  mixing  is  thought  to  be  an  important  (although 
perhaps  not  dominant)  process  occurring  frequently  in  the  Athabasca  Oil 
Sands  area.  During  daylight  hours  throughout  late  spring  to  early 
autumn,  thermal  mixing  will  also  be  Important;  however,  mechanical 
mixing  will  still  be  important  when  winds  at  plume  height  are  greater 
than  approximately  6 to  8 ms”^ . 

In  a sensitivity  and  validation  study  of  a Gaussian  frequency 
distribution  model,  Davison  et  al . (1981b)  showed  that  changing  Zq 
from  0.3  to  0.9  m caused  marked  changes  in  both  the  location  and 
magnitude  of  the  maximum  ground  level  concentration  (GLC)  values.  The 
changes  were  functions  of  the  source  character  1st i cs,  the  thermal 
stability,  and  the  wind  speed.  For  example,  at  a downwind  distance  of 
5 km  in  mechanically  dominated  mixing,  the  sector-averaged  GLC  values 
were  increased  by  over  40^  when  was  increased  from  0.3  to  0.9  m. 

This  was  a greater  effect  than  changing  the  wind  speed  from  10  to 
1 5 ms“^ . 

In  stable  conditions,  the  effect  of  changing  Zq  from  0.3  to 
0.9  m was  to  increase  the  GLC  value  by  more  than  a factor  of  two  and  to 
change  the  location  of  the  maximum  by  many  kilometres.  The  sensitivity 
study  also  showed  that  adopting  a value  for  Zq  of  0.9  m largely 
removed  any  systematic  discrepancies  between  predicted  and  observed  GLC 
values  within  the  limits  of  the  available  data. 
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3.  CHARACTERISTICS  OF  THE  MINISONDE  DATA  BASE 

3.1  DATA  SOURCES 

The  minisonde  data  base  used  to  determine  Zq  was  collected 
in  the  Athabasca  Oi I Sands  area  from  1975  to  1979.  Many  of  the 
releases  took  place  during  intensive  field  studies  in  1976  and  1977. 
Both  single  and  double  theodolite  readings  were  taken;  most  of  the 
single  theodolite  measurements  were  taken  prior  to  1977.  Data  from 
about  2000  minisondes  released  during  dayiight  hours  were  available. 

For  the  double  theodolite  releases,  observational  intervals  were  15  to 
30  s which  provided  spatial  resolutions  of  about  50  m. 

3.2  DESCRIPTION  OF  SITES 

The  Athabasca  Oil  Sands  area  Is  character i zed  by  a river 
valley  within  a region  of  rolling  terrain.  The  major  topographical 
features  are  the  Birch  Mountains  running  southwest  to  northeast  about 
40  km  northwest  of  the  site,  Stoney  Mountain  about  40  km  to  the  south, 
and  a gradual  rise  to  Muskeg  Mountain  in  the  east.  The  ground  cover  is 
a mixture  of  white  spruce  and  aspen  and  open  black  spruce  stands  within 
fen  and  bog  areas  (Thompson  et  al . 1978). 

The  minisonde  data  examined  in  this  report  were  from  two 
sites.  Syncrude  and  Lower  Syncrude  (see  Figure  1).  The  Syncrude 
minisonde  site  is  located  about  50  m southeast  of  the  Syncrude  plant. 

To  the  south  through  west  to  northeast  (170°  to  040°)  within  about 
1.5  km  from  the  release  site  are  plant  buildings  of  one  to  several 
stories  in  height  and  the  main  stack  which  is  over  180  m high.  Beyond 
about  2 km  are  scattered  strip  mines  free  of  vegetation  intermixed  with 
semi-open  white  spruce  and  aspen  forest.  Southeast  of  the  release 
site,  the  land  drops  away  slowly  with  a slope  of  about  10m  in  4 km  and 
is  covered  by  semi-open  white  spruce  and  aspen  forest  ranging  in  height 
up  to  10  m. 

The  Lower  Syncrude  site  is  located  within  the  Athabasca  River 
valley  flood  plain  about  200  m west  of  the  river  and  1 km  east  of  the 
west  bank  of  the  valley.  The  bank  rises  about  60  m in  200  m and  the 
the  bank  axis  is  oriented  north-northwest  to  south-southeast  for 
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Figure  1.  Map  of  Syncrude  and  Lower  Syncrude  mlnisonde  release  sites  in  late  1976.  Contour 
Intervals  are  10  m.  Horizontal  scale  Is  in  metres. 
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several  kilometres  in  either  direction  from  the  release  site.  To  the 
north  for  1 km  and  northwest  for  several  kilometres  is  a flat 
bog/fenland.  Similar  vegetation  exists  to  the  east  for  200  m;  beyond 
that  lies  the  river  which  is  approximately  800  m wide.  About  100  m to 
the  southeast  and  south  of  the  release  site  is  a smal 1 lake  about  300  m 
wide;  beyond  the  lake  is  a grouping  of  low  buildings  typically  10  m or 
less  in  height. 

3.3  UNCERTAINTIES  IN  THE  MIN  I SONDE  MEASUREMENTS 

Uncertainties  inherent  in  calculations  of  wind  speed  and 
direction,  temperature  lapse  rate,  and  balloon  height  from  double 
theodolite  minisonde  techniques  were  examined  by  Schaefer  and  Doswel I 
(1978)  and  extended  by  Netterville  and  Djurfors  (1979).  Measurement 
errors  were  shown  to  accumulate  with  time  since  release;  errors 
associated  with  readings  below  200  m (less  than  two  minutes  from 
release  with  typical  ascent  rates  of  2 ms"^ ) were  less  than  10^. 

Note  that  this  is  the  minimum  expected  error;  observer  error  and  lack 
of  instrument  resolution  will  increase  this  value. 

Atmospheric  turbulence  introduces  additional  uncertainties 
into  estimates  of  the  ensemble-average  profile  measurements.  Following 
Nettervi I le  and  Djufors  (1979), 

T 

= 1/T  (ay/u)2  y-R^Ct)  (1  - t/T)  dt  (1) 

o 

where  is  the  relative  wind  error,  T is  the  averaging  time  between 
consecutive  position  fixes,  <^y/U  is  the  along-wind  turbulent 
intensity,  and  R|_  is  the  Lagrangian  autocorre  I at  i on  coefficient.  The 
typical  averaging  times  for  minisondes  are  significantly  less  than  the 
Lagrangian  integral  time  scales  and  so  the  error  estimate  is 
approxi mate  I y 

e n.  0.7  a /U 
u u 


(2) 
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The  value  of  cTy/U  is  a function  of  roughness  length  and  stability 
and  also  involves  low-frequency  contributions  which  do  not  obey 
Mon i n-Obukhov  scaling  (Panofsky  1973).  A typical  value  in  strong  winds 
(neutral  conditions)  at  a height  of  100  m with  Zq  = 1 m is 


k/(u#  In  0.2 

where  k is  von  Karman’s  coefficient.  Combining  (2)  and  (3)  and  noting 
that  the  error  estimates  associated  with  the  minisonde  readings  are 
independent  leads  to  a probable  error  of  about  20^. 

Much  of  the  mini  sonde  data  produced  in  the  Athabasca  Oil 
Sands  area  has  been  examined  by  Davison  and  Leavitt  (1979)  who  cited 
several  cases  of  profiles  that  were  likely  incorrect  based  on  known 
meteorological  conditions  and  on  comparisons  with  other  profiles.  They 
also  estimated  errors  in  wind  speed  to  be  about  20^. 

Some  of  the  minisonde  data,  especially  prior  to  1977,  are 
based  upon  single  theodolite  measurements  with  an  assumed  rise  rate. 

The  error  associated  with  single  theodolite  measurements  over  the  first 
several  hundred  metres  is  probably  largely  dependent  upon  the  precision 
of  balloon  inflation  by  the  min  I sonde  technician. 
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4.  ROUGHNESS  LENGTH  CALCULATIONS 

4.1  PROFILE  ANALYSIS 

Mon i n-Obukhov  similarity  theory  has  proven  to  be  effective 
for  interpreting  atmospheric  boundary  layer  wind  profiles.  The  wind 
shear  can  be  expressed  as  follows  (Businger  1973): 

\ = kZ/u#  3U/aZ  (4) 

= (1  - 15  Z/L)"^/^  for  Z/L  < 0 (unstable)  ^ (5) 

= 1 +5  Z/L  for  Z/L  > 0 (stable)  (6) 

where  k is  von  Karman's  constant  and  L is  the  Mon i n-Obukhov  length. 
Explicit  expressions  for  the  wind  profile  (Paulson  1970)  are,  for 

Z/L  < 0: 

U/u.  = (In  Z/Zq  - H'l  )/k  (7) 

= In  [d+X)2(l+x2)/8]  -2  tan"'(X)  + H/2  (8) 

where 

X = = (1-15  Z/L)'/'‘ 

and  for  Z/L  > 0: 

U/u,  = (In  Z/Zq  - 1<2)/k 

where 

= -4.7  Z/L  (11) 

Because  fluxes  were  not  directly  measured,  the  bulk  Richardson  number 
was  used  to  estimate  Z/L  (following  Arya  1982) 


Z/L  = Ri 

when  Ri  < 0 

(12) 

Z/L  = Ri/(1-5  Ri  ) 

when  Ri  > 0 

(13) 

(9) 


(10) 
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where 

Ri  = g/0  (A0/AZ)  (14) 

and  0 is  potential  temperature,  Z is  the  geometric  mean  height,  and 
is  the  wind  speed  nearest  the  150  m level.  These  equations  have 
been  found  to  be  generally  valid  to  heights  as  high  as  150  to  200  m 
(Lumley  and  Panofsky  1964). 

The  Richardson  number  and  the  Mon i n-Obukhov  length  were 

calculated  from  the  minisonde  data.  The  Z was  estimated  by  a least 

o 

squares  fit  of  U versus  (in  Z -T)  for  each  profile.  Since  the  first 
temperature  value  was  typically  at  50  or  60  m (surface  data  were  not 
provided)  the  diabatic  influence  for  convective  conditions  could  be 
significantly  underestimated,  even  through  diabatic  effects  may  be 
large  enough  to  mask  Z^  effects  on  the  profiles.  Minisonde  selection 
criteria  and  prof i I e-by-prof i I e examination  were  used  to  ensure  that 
candidate  profiles  were  indeed  mechanically  dominated  and  coupled  to 
the  surface. 

4.2  MINI  SONDE  SELECTION  CRITERIA 

The  2000  minisonde  profiles  used  in  the  study  were  subjected 
to  a selection  process  to  produce  well-behaved  candidate  profiles  for 
roughness  length  calculations.  Required  profile  attributes  included: 

1.  Wind  speeds  in  access  of  about  5 ms  * and  temperatures 
decreasing  with  height  to  ensure  that  the  flow  was  coupled 
to  the  surface;  and 

2.  Wind  speeds  increasing  approximately  logarithmically  with 
height  after  compensation  for  diabatic  effects. 

Three  data  points  were  considered  a minimum  number  to  define 
a profile.  A constraint  was  also  required  on  the  wind  turning  with 
height  as  a large  turning  could  indicate  the  presence  of  a stable  layer 
and  hence  decoupling  of  the  winds  from  the  surface. 

The  data  selection  criteria  that  produced  these  requirements 

1.  At  least  three  measurements  below  a maximum  height; 

2.  Miriimum  wind  speeds  greater  than  a specified  value; 

3.  Monoto , .1  ca  I I y increasing  wind  speeds  and  decreasing 
temperature  with  height; 


i nc I uded : 
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4.  variability  constraint  on  wind  direction,  wind  speed,  and 
temperature  profiles; 

5.  temperature  lapse  rate; 

6.  date  of  data  collection  (to  stratify  by  season). 

Table  1 shows  the  number  of  qualifying  profiles  for  various 

combinations  of  selection  criteria.  The  requirements  of  a minimum  wind 
speed  of  5 ms"^  at  all  levels  with  reasonable  constraints  for  a 
uniform  profile  of  wind  speed,  wind  direction,  and  temperature  resulted 
In  a very  large  reduction  In  the  number  of  qualifying  profiles.  This 
reduction  Is  consistent  with  the  c I I mato I og I ca I I y light  wind  speeds  of 
the  area  (Long  ley  and  Janz  1978).  Further  corrstraints  on  the 
temperature  lapse  rate  and  the  maximum  height  for  the  lowest  three 
measurements  were  considered  to  be  prudent  for  estimating 
roughness  lengths.  The  selection  criteria  adopted  for  the  roughness 
length  calculations  are  shown  In  Table  1 and  consisted  of  150  m maximum 
height  and  a lapse  rate  within  0.5°C/100  m of  adiabatic.  Note  that 
lapse  rate  had  little  effect  on  the  number  of  profiles  chosen.  The 
criteria  which  most  reduced  the  data  set  were  5 ms”"'  minimum  wind 
speed  and  the  three-level  profile.  No  attempt  was  made  to  stratify  by 
season  because  of  the  small  number  of  selected  profiles.  While  the 
small  number  of  profiles  might  have  increased  the  uncertainty  of  the 
mean  Zq  values,  rigorous  selection  ensured  that  only  well  behaved 
prof i I es  were  used. 

4.3  ROUGHNESS  LENGTH  ESTIMATES 

Roughness  length  estimates  for  profiles  meeting  the 
selection  criteria  are  shown  in  Table  2.  Included  are  estimates  with 
no  diabatic  effects,  with  diabatic  effects  but  no  displacement  heights, 
and  with  diabatic  effects  for  a range  of  assumed  displacement  heights. 

Bulk  Richardson  numbers  near  zero  in  Table  2 resulted  In 
values  corrected  for  diabatic  effects  being  approximately  equal  to 
neutral  Zq  values.  When  gradient  Richardson  numbers  were  used  for 
comparison,  as  in  Ming  et  al.  (1983),  increased  variability  in  both 
R|  and  Zq  were  noted.  Gradient  Rj  tended  to  be  larger  in 
absolute  value,  with  some  exceeding  the  assumed  critical  value  of  0.20. 
Thus,  the  bulk  formulation  for  Rj  was  used  exclusively. 
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Table  1.  The  number  of  qualifying  minisonde  profiles  from  a total  of 
2000  as  a function  of  selection  criteria.^ 


Min.  Height 
(m) 

Lapse  Rate 
(°C/100m) 

Date 

Qua  1 i f y i ng  Prof i 1 es 

200 

-1  + 0.5 

Al  1 

33 

200 

-1  + 0.5 

1977 

18 

150 

<0 

Al  1 

37 

150 

-1  + 0.5 

Al  1 

20  b 

150 

-1  + 0.5 

1977 

6 

150 

-1  + 0.25 

Al  1 

15 

150 

-1  + 0.1 

Al  1 

9 

150 

-1  + 0.1 

1977 

3 

^ Additional  criteria  applied  to  all  subsets  were: 


( i ) 

( i i ) 
( i i i ) 
( i v) 

(v) 


wind  speed  at  lowest  level  (about  50  m)  greater  than 
5 ms"^ ; 

wind  speeds  monoton i ca I I y increasing  with  height; 
wind  direction  variation  less  than  15°; 
temperature  within  1 °C  of  the  temperature  derived  for 
that  level  from  the  mean  linear  lapse  rate;  and 
wind  speed  within  ^ 1 ms”^  of  the  wind  speed  derived 
for  that  level  from  the  mean  logarithmic  profile. 
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Finally  selected  criteria. 


Zp,  with  dlabatlc  correction  (m) 

Profile  R|^  Wind  Zq  No 

XlO-2  (m)  Direction  Neutral  Displacement  DISI^  DIS2C  DIS1/2C 
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DISl  refers  to  the  application  of  physically  reasonable  displacement  heights  as  a function  of 
direction.  The  assumed  displacement  heights  for  Lower  Syncrude  are  5 m for  the  wind  direction  range 
(150®,  340®)  and  zero  elsewhere.  The  assumed  displacement  heights  for  the  Syncrude  site  are  5 m for 
(180®,  250°),  10  m for  (250®,  20°),  and  zero  elsewhere. 

DIS2  and  DISl/2  refer  to  times  2 and  times  1/2  those  assumed  displacement  heights. 
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While  diabatic  effects  were  very  small,  terrain  effects  on 
minisonde  profiles  were  larger.  At  Lower  Syncrude,  profile  6 could 
have  been  influenced  by  the  valley  wall  located  about  1 km  westward, 
especially  considering  the  apparently  stable  boundary  layer.  It 
appears  to  be  the  only  profile  at  Lower  Syncrude  for  which  displacement 
heights  were  likely  significant. 

At  the  Syncrude  site,  winds  blowing  over  the  plant  could  have 
been  affected  by  the  increased  roughness  and  might  have  experienced  a 
displacement  height.  Physically  reasonable  displacement  heights 
reduced  the  Zq  estimates  by  an  average  of  about  15^.  Calculating 
estimates  of  the  displacement  heights  from  the  data  was  considered 
Inappropriate  based  upon  the  probable  error  of  the  estimates  as 
discussed  in  Section  5. 

A summary  of  average  Zq  estimates  by  release  location  Is 
presented  in  Table  3.  Because  of  large  variations  in  Zq  (several 
orders  of  magnitude),  both  arithmetic  and  logarithmic  averages  are 
presented.  Arithmetic  averages  of  Zq  were  near  1 m at  Lower  Syncrude 
and  near  8 m at  Syncrude  with  small  differences  among  estimates  using 
neutral,  diabatic  and  displacement  height  effects.  Arithmetic 
differences  In  Zq  between  the  two  sites  were  large  and  were  of  the 
same  order  as  the  standard  deviations.  Site-to-site  variations  In  Zq 
were  also  large  when  logarithmic  averages  were  used.  It  Is  evident 
that  Zq  values  from  the  selected  profiles  at  Syncrude  had  less 
variability  than  at  Lower  Syncrude,  as  indicated  by  the  relatively 
small  difference  between  arithmetic  and  logarithmic  averages  at  the 
Syncrude  site  and  by  the  ratio  of  standard  deviations  to  arithmetic 
averages.  Logarithmic  averages  of  Zq  at  Syncrude  were  near  4 m.  It 
should  be  noted  that  the  use  of  standard  deviations  does  not  Imply 
normally  distributed  Zq  values  at  either  site;  rather,  the  standard 
deviations  give  an  indication  of  the  variability  of  the  data. 

Spatial  differences  In  Zq  between  the  two  sites  were 
evaluated  using  the  non-parametr i c Mann-WhItney  U test.  This  test 
makes  no  assumptions  about  the  distribution  of  the  Zq  values,  ranking 
the  values  from  highest  to  lowest.  Very  large  or  very  small  values  of 
the  statistic  U Imply  a separation  of  the  ordered  values  and  Indicate  a 


Table  3.  Average  roughness  length  estimates  by  release  location. 
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difference  between  the  population  distributions.  Using  sample  sizes  of 
6 and  11  for  Lower  Syncrude  and  Syncrude,  U values  of  6 and  60  were 
found.  Using  the  lower  of  the  values,  U = 6,  as  the  test  statistic 
(Seigel  1956)  the  one-tailed  Mann-Whitney  U test  indicates  that  is 
smaller  at  Lower  Syncrude  than  at  Syncrude,  at  the  ]%  level. 

Wind  direction  dependencies  were  also  investigated,  but  no 
differences  were  found.  A I arger  data  set  is  likely  required  to 
resolve  this. 
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5.  DISCUSSION  OF  THE  ROUGHNESS  LENGTH  ESTIMATES 

The  surprisingly  large  values  of  roughness  length  presented 
in  Section  4 are  discussed  below  in  terms  of  uncertainties  in  the 
estimates,  comparisons  with  other  sites,  independent  support  from  other 
measurements  in  the  oil  sands  area,  and  a possible  stability 
dependence. 

5.1  UNCERTAINTIES  IN  THE  ESTIMATES 

The  most  probable  error  ( 5F)  of  a derived  parameter  (F) 
which  is  a function  of  Xj  constituent  measurements  Is  given  by  (see, 
for  example,  Baird  1962): 

(5F  = (6X.  3F/3X.  )^l 

I V=1  ' ' > f>5) 

Applying  this  formal  methodology  to  the  roughness  length  estimates  is 
complicated  by  the  least  squares  fitting  of  profiles  and  by  the  highly 
non-linear  effect  of  velocity  perturbations  on  the  Zq  estimates. 

An  alternative  procedure  presented  by  Blanc  (1983)  involves 
explicit  calculation  of  the  parameter  F for  the  error  estimates  of 
Xj.  His  procedure  Is  especially  convenient  for  computer-based 
analysis  since  the  error  perturbations  can  be  treated  by  the  same  code 
used  for  the  actual  data  analysis.  Blanc's  perturbation  approach  Is 
given  by: 


6F 


xi 


6F  = 


f|+  |Fxi  - F 


1/2 


(16) 


(17) 


where  F^j  is  the  value  of  F calculated  for  a positive  perturbation 
of  the  Xj  constituent  measurement.  For  this  application,  F=  Z^. 
Constituent  measurements  were  wind  speed,  temperature  gradient,  and 
ba I loon  he i ght. 
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A "base  case"  was  chosen  and  both  positive  and  negative 
perturbations  were  applied.  The  base  case  was  defined  as  having 
profile  measurements  at  50,  100,  and  150  m with  U(50)  = 5 ms“^  and 
a neutral  lapse  rate.  The  uncertainty  estimates  were  calculated  for 
Zq  values  of  1,  5,  and  10  m (from  which  the  wind  speeds  at  higher 
elevations  were  computed).  Wind  speed  uncertainties  were  20$  and 
temperature  and  balloon  height  uncertainties  were  taken  as  10$  (from 
Nettervi I le  and  Djurfors  1979).  Table  4 presents  a summary  of  results 
for  Zq  = 5 m and  a range  of  wind  speed  perturbat ions;  Table  5 
presents  uncertainties  for  a range  of  Zq.  Uncertainty  estimates 
detailing  the  contributions  to  the  total  probable  error  of  Individual 
measurements  (wind  speed,  height,  and  temperature  lapse  rate)  are  given 
in  the  appendix.  These  estimates  Indicated  that  errors  in  wind  speed 
were  dominant  for  wind  speed  perturbat ions  larger  than  about  5$. 

For  wind  speed  perturbat i ons  of  5 to  30$,  the  total  probable 
error  in  Zq  (from  Table  4)  is  shown  to  range  from  4 to  24  m,  for  an 
initial  Zq  of  5 m.  A factor  of  six  variation  in  the  perturbation 
results  in  a factor  of  six  variation  in  the  probable  Zq  error.  The 
error  in  Zq  is  sensitive  to  wind  speed  perturbat ion . 

For  wind  speed  perturbat ions  of  20$  and  initial  Zq  ranging 
from  1 to  10  m (Table  5),  the  range  in  Zq  error  was  much  smaller.  In 
fact,  in  the  Zq  range  from  1 to  10  m,  the  error  was  nearly  constant 
at  1 1 m.  This  was  expected  since  the  wind  speed  component  error 
appears  to  dominate  and  the  wind  speed  perturbation  was  held  constant. 
Note,  however,  that  the  error  in  In  Zq  did  not  remain  constant  but 
decreased  with  increasing  In  Zq. 

Roughness  length  estimates  in  Tables  2 and  3 can  be  compared 
to  error  estimates  In  Tables  4 and  5.  Standard  deviations  of  Zq 
(from  Table  3)  ranged  from  4 to  7 m with  most  near  7 m.  Probable 
errors  in  Zq  (from  Table  4)  ranged  from  3.7  m (5$  wind  speed 
perturbat ion ) to  24  m (30$  wind  speed  perturbation);  a standard 
deviation  of  7 m corresponded  to  a perturbation  of  about  13$.  That 
this  value  (13$)  is  less  than  the  20$  error  suggested  by  Davison  and 
Leavitt  (1979)  and  by  the  a^^/U  analysis  suggests  that  careful 
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Table  4.  Roughness  length  uncertainty  estimates  as  a function  of 

perturbation  in  U(Z)  for  neutral  profiles  with  U(50)  = 5 ms"^ 
and  Zq  = 5 m.  Total  error  is  the  total  probable  error  assuming 
independence  of  constituent  errors.  Uncertainty  estimates  are 
differences  in  metres  from  Zq  = 5 m (for  Zq). 


Perturbed 

Zo 

Perturbat ion 

i%) 

lnZ„ 

Perturbation  {%) 

Parameter 

5 

10 

20 

30 

5 

10 

20 

30 

U(50) 

1 .2 

3.5 

6.7 

9.4 

0.36 

0.73 

1 .5 

2.0 

unoo) 

0.35 

1 .4 

6.2 

14 

0.066 

0.25 

0.80 

1.3 

U(150) 

1 .7 

3.2 

4.7 

16 

0.33 

0.60 

0.55 

1 .4 

Errors  due  to 
Z and  dT/dZ 
(10^  pertur- 
bation) 

2.7 

2.7 

2.7 

2.7 

0.60 

0.60 

0.60 

0.60 

Total  error 
i n ZQ(m) 
or  IhZq 

3.7 

5.7 

11 

24 

0.78 

1.1 

1 .9 

2.9 
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Table  5.  Roughness  length  uncertainty  estimates  as  a function  of  Zq 

for  neutral  profiles  with  U(50)  = 5 ms“^  and  20^  perturbation 
in  U(Z).  Total  error  is  the  total  probable  error  assuming 
independence  of  constituent  errors.  Uncertainty  estmates  are 
differences  In  meteres  (for  Zq)  from  Indicated  Zq  values). 


Perturbed 

^o 

Z^(m) 

In  Zo 
Z„(m) 

Parameter 

1 

5 

10 

1 

5 

10 

U(50) 

4.3 

6.7 

7.5 

2.2 

1 .5 

0.83 

U(IOO) 

7.9 

6.2 

5.0 

2.1 

0.80 

0.41 

U(150) 

9.5 

4.7 

6.1 

2.3 

0.55 

0.51 

Errors  due  to 
Z and  dT/dZ 
( 10^  perturbat ion ) 

1 .3 

2.7 

3.2 

5.4 

0.60 

0.31 

Total  error  in 

13 

11 

1 1 

6.6 

1 .9 

1.1 

ZQ(m)  or  InZQ 
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selection  of  candidate  profiles  helps  to  reduce  apparent  errors  in 
measurement.  The  scatter  in  Zq  estimates  in  Table  3 can  largely  be 
accounted  for  by  reasonable  uncertainties  in  profile  measurements. 

5.2  COMPARISON  OF  Zq  VALUES  MEASURED  AT  OTHER  SITES 

The  Zq  estimates  from  this  study  can  be  compared  to  values 
In  other  regions.  For  example,  Korrell  et  al.  (1982)  analysed  profile 
data  between  10  and  50  m from  the  Boulder  300-m  instrumented  tower. 

This  site  is  reasonably  flat  with  terrain  slopes  near  2%  and  low 
vegetation  with  occasional  trees  and  houses  within  a radius  of  about 
3 km.  Zq  was  found  to  be  di rect ion-dependent  with  values  ranging 
from  4 to  35  cm.  Ming  et  al.  (1983)  analysed  profiles  at  three  towers 
in  New  England.  At  one  site,  surrounded  by  tall  and  Irregular  trees  in 
rolling  terrain,  Zq  values  ranged  from  1 to  1 1 m based  on  data  at  46 
and  99  m.  At  a second  site.  In  forest  and  rolling  farmland,  Zq 
values  were  near  1 m (21,  46,  and  99  m measurements),  while  at  a third 
site,  with  water,  woods,  farmland  and  buildings  in  different 
directions,  Zq  varied  from  less  than  1 cm  to  about  250  cm  (10,  43, 
and  114  m measurements). 

5.3  THE  IMPLICATION  OF  DISSIPATION  MEASUREMENTS  ON 
ROUGHNESS  LENGTH 

Independent  supporting  evidence  for  the  existence  of 
relatively  large  Zq  can  be  found  in  height-dependent  turbulent 
dissipation  values  measured  in  the  Athabasca  Oil  Sands  area. 

Typically,  the  dimensionless  dissipation  rate  ej  decreases  rapidly 
with  height  near  the  surface  and  then  is  approximately  constant  within 
the  convectively  mixed  layer  (see,  for  example,  Kaimal  et  al . 1976); 

The  dissipation  rate  is  given  by: 

"^D  = "^T/(g  Qq)  (18) 

where 

e = turbulent  energy  dissipation; 

T = mean  temperature;  and 
Qq  = surface  kinematic  heat  flux. 
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A constant  requires  e to  increase  with  height  as  T decreases.  In 

the  oil  sands  area,  however,  Davison  and  Grandia  (1979)  found  e near 
plume  level  on  clear  sunny  days  to  consistently  decrease  with  height. 
They  concluded,  therefore,  that  mechanical  mixing  was  Important  at 
plume  heights  even  In  the  presence  of  significant  surface  heat  fluxes. 
Venkatram  (1980),  on  the  other  hand,  suggests  that  large  surface  heat 
fluxes  In  the  oil  sands  area  Imply  that  free  convection  should  occur 
frequently.  These  two  pieces  of  evidence  can  be  reconciled  If  Is 
large.  In  that  case,  ep  would  decrease  with  height  as  surface 
effects  are  observed  at  higher  levels,  and  therefore  free  convection 
would  not  be  expected  to  occur  frequently  even  with  large  heat  fluxes. 

5.4  A POSSIBLE  STABILITY  DEPENDENCE  FOR  ROUGHNESS  LENGTH 

The  application  of  (large)  Zq  values  from  this  study  to  air 
quality  modelling  studies  In  the  Athabasca  Oil  Sands  area  requires 
clarification.  In  slightly  unstable  conditions  the  data  Indicate  Zq 
values  ranging  from  1 to  8 m.  These  values  seem  reasonable  when  the 
contribution  to  Zq  from  form  drag  Is  Important.  In  this  case  Zq 
will  be  representat I ve  of  conditions  over  a wide  area  (up  to  a 
kilometre  or  more)  Including  the  effects  of  small  terrain  features. 

The  application  of  a large  Z^  Is  more  uncertain  In  stable 
conditions  for  several  reasons.  First,  air  tends  to  flow  around 
obstacles  rather  than  over  them,  causing  a decoupling  of  the  air  from 
the  surface  and  therefore  a situation  In  which  elevated  winds  are  not 
determined  by  the  underlying  terrain.  Second,  air  quality  model 
sensitivity  studies  In  the  area  by  Davison  et  al . (1981b)  suggest  Z^ 
to  be  of  the  order  of  0.5  to  1 m,  neglecting  dependence  on  stability. 
The  use  of  Zq  8 m would  Introduce  minor  changes  In  their  estimated 
GLC  values  In  unstable  conditions  but  would  Impose  a large,  systematic 
bias  to  large  GLC  values  In  stable  conditions.  Finally,  the  number  of 
stable  profiles  analysed  In  this  study  was  too  small  and  the  results 
were  too  variable  to  suggest  a new  Zq  estimate  In  stable  conditions. 
The  net  result  In  stable  conditions  Is  suggested  to  be  a smaller  value 
of  Zq  since  the  effects  of  form  drag  are  reduced. 
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The  results  of  this  study  suggest  that,  when  form  drag  is  an 
important  determinant  of  Z^,  the  roughness  length  may  be  a function 
of  stability.  This  could  occur  because,  while  the  ground  cover 
component  of  is  expected  to  vary  little  with  stability,  the 
effects  of  form  drag  might  be  expected  to  vary  substantially. 
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6.  CONCLUSIONS  AND  RECOMMENDATIONS 

6.1  CONCLUSIONS 

Roughness  lengths  were  calculated  from  profiles  of  mini  sondes 
released  in  the  Athabasca  Oil  Sands  area  from  1975  to  1979.  A rigorous 
selection  process  to  obtain  well-behaved  wind  and  temperature  profiles 
resulted  in  only  20  suitable  profiles  at  two  sites  from  an  original 
database  of  about  2000  minisonde  releases. 

Based  on  measurement  levels  typically  near  50,  100, and  150  m, 
mean  values  of  were  found  to  range  from  1 m at  Lower  Syncrude  to 
8 m at  the  Syncrude  plant  site.  The  Mann-Whitney  U test  showed  these 
differences  to  be  significant  at  the  \%  level.  These  values  were 

comparable  to  other  values  in  similar  terrain  quoted  in  the  literature 
and  were  consistent  with  height-dependent  dissipation  rates  measured  in 
the  area.  No  differences  with  wind  direction  were  found.  It  was 
suggested  that  Z^  may  be  stability-dependent  in  terrain  where  form 
drag  is  important,  such  that,  in  stable  conditions,  the  effective 
roughness  length  is  much  smaller  due  to  the  decreased  effectiveness  of 
form  drag. 

An  error  analysis  showed  that  the  observed  variations  in  the 
Zq  estimates  were  similar  to  the  estimated  uncertainties  in  the 
constituent  measurements.  It  was  also  shown  that  uncertainties  in  the 
wind  values  were  much  nnor  e important  for  Zq  calculations  than 
uncertainties  in  the  balloon  height  or  in  the  temperature  gradient. 

6.2  RECOMMENDATIONS 

1.  It  is  recommended  that  Zq  be  determined  from  profiles 
from  the  150  m tower  located  at  Lower  Syncrude.  The 
tower  profiles  should  provide  a larger  data  base  than  the 
minisondes  and  should  have  smaller  associated  measurement 
errors  (especially  wind  speed).  The  tower  data  analysis 
would  emphasize  Zq  calculation  in  slightly  unstable 
conditions  and  should  provide  Zq  estimates 
representat i ve  of  conditions  in  the  Athabasca  River 
valley. 
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2.  It  is  recommended  that  a catalog  of  Zq  values  be 
compiled  for  Alberta.  This  task  would  include 
synthesizing  Zq  values  where  they  exist  and  calculating 
them  from  existing  profile  data  in  regions  where  tower  or 
other  atmospheric  sounding  data  are  available.  An 
instrumented  research  aircraft  might  be  used  in  areas 
where  no  profile  data  exist.  The  catalog  Is  recommended 
because  of  the  spatial  and  stability  dependence  of  Zq 
suggested  by  the  results  of  this  study. 
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8.  APPENDIX 

8.1  ROUGHNESS  LENGTH  UNCERTAINTY  ESTI^TES 

This  section  documents  the  contributions  to  total  probable 
error  of  the  individual  constituent  profile  measurements.  Probable 
errors  for  Zq  (differences  from  the  original  value)  are  in  metres. 

The  overbar  indicates  the  mean  of  a positive  and  negative  perturbation 
for  a single  parameter.  The  total  probable  error  in  Zq  and  In  Zq 
is  given  at  the  bottom  of  each  table. 
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Table  6.  Roughness  length  uncertainty  estimates  for  Zq  = 5 m, 
U(50)  = 5 ms“^,  and  a neutral  profile. 


Perturbed 

Perturbation 

Parameter 

% 

«Zoi 

«Zoi 

6 I hZq  j 

5 I n Zq  j 

U (50) 

“20 

8.9 

1.1 

U (50) 

+20 

4.4 

6.7 

2.1 

1 .5 

U (100) 

-20 

6.7 

0.86 

U (100) 

+20 

5.6 

6.2 

0.74 

0.80 

U (150) 

-20 

0.39 

0.081 

U (150) 

+20 

9.1 

4.7 

1 .0 

0.55 

Z (50) 

-10 

1 .4 

0.33 

Z (50) 

+ 10 

1 .7 

1 .6 

0.30 

0.31 

Z (100) 

-10 

0.068 

0.013 

Z (100) 

+ 10 

0.058 

0.063 

0.012 

0.012 

Z (150) 

+ 10 

1 .2 

0.21 

Z (150) 

+ 10 

0.86 

1 .0 

0.19 

0.20 

dT/dZ 

-10 

2.6 

0.72 

dT/dZ 

+ 10 

1 .3 

1 .9 

0.22 

0.47 

«Zo  = 11 
s InZo  = 1 .9 
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Table  7.  Roughness  length  uncertainty  estimates  for  Zq  = 1 m, 
U(50)  = 5 ms”^ , and  a neutral  profile. 


Perturbed 

Perturbation 

Parameter 

% 

^Zoi 

SInZoi 

<5 1 n Zq  j 

U (50) 

-20 

7.7 

2.1 

U (50) 

+20 

0.96 

4.3 

2.3 

2.2 

U (100) 

-20 

9.2 

2.3 

U (100) 

+20 

6.5 

7.9 

2.0 

2.1 

U (150) 

-20 

12 

2.5 

U (150) 

+20 

7.5 

9.5 

2.1 

2.3 

Z (50) 

-10 

0.41 

0.48 

Z (50) 

+ 10 

0.58 

0.49 

0.44 

0.46 

Z (100) 

-10 

0.048 

0.044 

Z (100) 

+ 10 

0.038 

0.043 

0.037 

0.040 

Z (150) 

-10 

0.41 

0.33 

Z (150) 

+ 10 

0.26 

0.34 

0.28 

0.31 

dT/dZ 

-10 

1.1 

10.0 

dT/dZ 

+ 10 

1 .2 

1 .1 

0.75 

5.4 

(SZq  = 13.0 

6|nZo  = 


6.6 
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Table  8.  Roughness  length  uncertainty  estimates  for  Zq  = 10  m, 
U(50)  = 5 ms”^,  and  a neutral  profile. 


Perturbed 

Perturbat ion 

Parameter 

% 

<^^oi 

<5^0  i 

6lnZoi 

5 1 nZo i 

U (50) 

-20 

8.6 

0.62 

U (50) 

+20 

6.5 

7.5 

1 .0 

0.83 

U (100) 

-20 

5.7 

0.45 

U (100) 

+20 

4.4 

5.0 

0.36 

0.41 

U (150) 

-20 

3.2 

0.38 

U (150) 

+20 

9.0 

6.1 

0.64 

0.51 

Z (50) 

-10 

2.3 

0.26 

Z (50) 

+ 10 

2.7 

2.5 

0.23 

0.25 

Z (100) 

-10 

0.014 

0.0014 

Z (100) 

+ 10 

0.010 

0.012 

0.0010 

0.0012 

Z (150) 

-10 

1 .6 

0.15 

Z (150) 

+ 10 

1 .3 

1 .5 

0.14 

0.14 

dT/dZ 

-10 

1 .6 

0.18 

dT/dZ 

+ 10 

0.87 

1 .3 

0.083 

0.13 

fiZo  = 11 

6 InZo  = 


1.1 
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Table  9.  Roughness  length  uncertainty  estimates  for  = 5 m, 
U(50)  = 5 ms"^,  and  a neutral  profile. 


Perturbed 

Perturbat ion 

Parameter 

% 

«Zoi 

«Zoi 

SinZoj 

6lnZoi 

U (50) 

-5 

1.9 

0.33 

U (50) 

+5 

1 .6 

1 .8 

0.40 

0.36 

U (100) 

-5 

0.27 

0.053 

U (100) 

+5 

0.42 

0.35 

0.080 

0.066 

U (150) 

-5 

1.5 

0.34 

U (150) 

+5 

1 .9 

1 .7 

0.32 

0.33 

Z (50) 

-10 

1 .4 

0.33 

Z (50) 

+ 10 

1 .7 

1 .6 

0.30 

0.31 

Z (100) 

-10 

0.068 

0.013 

Z (100) 

+ 10 

0.058 

0.063 

0.012 

0.012 

Z (150) 

-10 

1.2 

0.21 

Z (150) 

+ 10 

0.86 

1 .0 

0.19 

0.20 

dT/dZ 

-10 

2.6 

0.72 

dT/dZ 

+ 10 

1 .3 

1 .9 

0.22 

0.47 

6Zo  = 3.7 

S lnZ„  = 


0.78 
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Table  10.  Roughness  length  uncertainty  estimates  for  Zq 
U(50)  = 5 ms“^ , and  a neutral  profile. 


5 m, 


Perturbed 

Perturbat ion 

Parameter 

% 

«Zoi 

«Zoi 

6 1 hZq j 

6 1 n Zq  j 

U (50) 

-10 

4.1 

0.60 

U (50) 

+ 10 

2.9 

3.5 

0.87 

0.73 

U (100) 

-10 

1 .4 

0.24 

U (100) 

+ 10 

1 .5 

1 .4 

0.26 

0.25 

U (150) 

-10 

2.3 

0.60 

U (150) 

+ 10 

4.2 

3.2 

0.60 

0.60 

Z (50) 

-10 

1 .4 

0.33 

Z (50) 

+ 10 

1 .7 

1 .6 

0.30 

0.31 

Z (100) 

-10 

0.068 

0.013 

Z (100) 

+ 10 

0.058 

0.063 

0.012 

0.012 

Z (150) 

-10 

1.2 

0.21 

Z (150) 

+ 10 

0.86 

1 .0 

0.19 

0.20 

dT/dZ 

-10 

2.6 

0.72 

dT/dZ 

+ 10 

1 .3 

1 .9 

0.22 

0.47 

6 Zq  = 5.7 

6 InZo  = 1.1 
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Table  11. 


Roughness  length  uncertainty  estimates  for  Zq  = 5 m, 
U(50)  = 5 ms”^ , and  a neutral  profile. 


Perturbed 

Perturbat ion 

— 

Parameter 

% 

«2oi 

«Zoi 

SinZoi 

«lnZoi 

U (50) 

-30 

14.0 

1 .3 

U (50) 

+30 

4.7 

9.4 

2.7 

2.0 

U (100) 

-30 

17.0 

1.5 

U (100) 

+30 

12.0 

14.0 

1 .2 

1 .3 

U (150) 

-30 

19.0 

1 .6 

U (150) 

+30 

14.0 

16.0 

1 .3 

1 .4 

Z (50) 

-10 

1 .4 

0.33 

Z (50) 

+ 10 

1 .7 

1 .6 

0.30 

0.31 

Z (100) 

-10 

0.068 

0.013 

Z (100) 

+ 10 

0.058 

0.063 

0.012 

0.012 

Z (150) 

-10 

1 .2 

0.21 

Z (150) 

+ 10 

0.86 

1 .0 

0.19 

0.20 

dT/dZ 

-10 

2.6 

0.72 

dT/dZ 

+ 10 

1 .3 

1 .9 

0.22 

0.47 

= 24 
6 InZo  = 


2.9 
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